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Abstract 

Coupled DGLAP equations involving singlet quark and gluon distributions 
are explored by a Taylor expansion at small x as two first order partial differ- 
ential equations in two variables : Bjorken x and t {t = log^).Th.e system of 
equations are then reduced to canonical form and the resultant equations are 
solved by the method of characteristics Analytic and semi-analytic solutions 
thus obtained are compared with the exact results and the range of validity 
obtained. 

PACS Nos: 12.38.-t;12.38.B_x;13.60.-r;13.60.Hb 



1 Introduction 



DGLAP equations [H |2l El El have been playing very important role in the study 
of nucleon structure function in deep inelastic regime. In the standard procedure, 
the X profile of the parton distribution at some low scales Qq are parametrized by 
comparing with data while their values at a higher scale are obtained by evolution 
of these equations, which are solved numerically. In the small x regime however, 
several analytical solutions^ O 13 El El of the DGLAP equations are available in 
the literature, which are in good agreement with HERA data [1011111. This suggests 
utility of such approaches in understanding the dynamics of evolution of quarks and 
gluons at small x. In recent years the present authors and their collaborators |12| 
ir^l HH US] have also been pursuing an approximate method of solution of DGLAP 
equations at small x with considerable phenomenological success 

In this paper we present solutions of the coupled singlet evolution equations 
applying the method of characteristics. By a Taylor series expansion of the gluon and 
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the singlet distribution that occur under the integral sign of the coupled equations, 
we convert the coupled integro-differential equations into two first order partial 
differential equations in two variables : Bjorken x and t{t = ln{Q'^/h?)). The 
resulting system of equations which are still coupled, is then reduced to canonical 
form by the introduction of a suitable new function. After reduction to the canonical 
forms, the equations can be solved by applying the method of characteristics [il^ IT7| 
ITH] under certain approximations valid to be at small x. 

In ^ we describe the formahsm, ^ is devoted to testing our prediction with 
the exact solution [12] and data [201111] , while in |2]we give our conclusion. 



2 Formalism 



2.1 Singlet coupled DGLAP equations in Taylor approxi- 
mated form 

The coupled DGLAP equations for quark singlet (S(a;,Q^)) and gluon {g{x,Q'^)) 
densities are P |2l 13 HI 



dlogQ^ I 9{x,Q'') 



asm_( P,, P,A ( j:{x,q 

2n I Pg, Pgg [ 9{x,Q'- 



where as (Q^) is the strong coupling constant, Pijs are the Altarelli-Parisi splitting 
functions and the symbol C?> stands for the usual Mellin convolution in the first 
variable already defined as 



a{x) (g) f{x) 



^ y \y, 



(2) 



Introducing the variable t = In^ and using the explicit forms of the splitting 
functions Pjj in LO, Eq.(0) can be written as 



dFi{x,t) Af 
dt t 



1 dz 



{3 + Aln{l-x)}Fi{x,t) + 2j^ - 



-2Fi{x, t)} + Inj !l dz {z" + (1 - zf) G (f , t 



X 



(3) 



and 



dG{x,t) Af 



dt 



+ dz 

J X 



l-Z 



+ (z{l-z) + 



t 

zG[%t)-G{x,t) f ^ \-z\^ix 



z 



+^ /' dz 

O Jx 



1 + (1 - ^) 



0. (4) 
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Here Af = Po = 11 — |n/ and as{t) = ^ . G{x,t){= xg{x,t)) is the gluon 
momentum distribution and F2{x,t) is the singlet structure function of the proton 
defined as 

"/ 

Fi{x, t) = xS(x, t) = xJ2 t) + U^, t)) , (5) 

2=1 

where qi{x,t) is the quark distribution of the i th flavour inside the proton. We 
write Eqs.© and (jlj) as 



dFi{x,t) Af 
dt t 



{3 + Aln{l - x)}F^{x, t) + It{x, t) + If (x, t) 







and 



dG{x,t) QAj 



dt 



li _ ^ + Inil - x) \ t) + F^ix, t) + t) 



(6) 



0, (7) 



where 



n{x,t) = 2j^^ 



z 



{l + z'')Fi{-,t)-2F,'{x,t) 



,x 



X 



mx,t) = -nf / dz[z' + il-z)']G{-,t) 



mx,t) 



2 /-i 



dz 



X 



and 



Iiix,t) 



1 (zG(^^,t)-G{x,t) 
dz ■' 



1 - z 



+ {z{l-z) + 



1- z 



G{^,t 



(9) 
(10) 

(11) 



To carry out the integrations in Eqs.ip irUI) . we introduce the variable u defined as 



X 



u = 1 — z and expand the argument — in F2*(-,t) and G{-,t) as a series. 



X 



X 



z 1 — u 



X^U^ = X + x^u^. 



k=0 k=l 

Using Eq.f IT^ we expand F2*(-,t) and G{-,t) in Taylor series as: 



Fi{-,t) = Fi{x,t)+xJ2u' 



k=l 



and 



G{-,t) =G{x,t)+xJ2 



dFjix^t) ^ ixEk=iu'yd'Fi{x,t) 
dx 2! dx"^ 

^dGix,t) (a;Er=iw')'^'G(x,t) 



u 



k=l 



dx 



+ 



2! 



dx'^ 



+ 



(12) 



(13) 



(14) 



The series Eqs. (II 3j) and (fT^ are convergent fH^ IT^ and hence at small x , we can 
approximate these by 



Fi{-,t)^Fi{x,t)+xJ2u' 
^ k=i 



■ dFiix,t) 
dx 



(15) 



and 

t) ^ G{x, t) + X V u'^^i^. (16) 

Using the above two Eqs. ()15|1 and ()16|) we carry out the integrations in z in Eqs. (jHJ- 
Uni). Neglecting terms O(x^) which is justified at small x , we get 



and 



mx,t) ^ {2x-3)Fi{x,t)+ x + 2xln- , (17) 

V xj ox 

If(x, t)^nf(l- -x] G{x, t) - !^ fsx - 3a; In -) ^^i^, (is) 
\ 2 J 2 \ xJ ox 

I2M ^ (-X + ^xlni) Fi{x,t) + Ix^^^i^, (19) 

II(x,t)^R'^{x)Gix,t) + P,ix)^^^, (20) 



where the function Pg{x) is 



P,(a;) = x(lni-H), (21) 



and R'g{x) is given by 

i?;(x) = (/ni-ii) + 2x. (22) 

Using Eqs.f imOUj) . we recast Eq.® and Eq.© as two first order partial differential 
equations in x and t in standard form: 

dFi{x,t) , 9G(M) a?f(M) >9G'(x,t) 

On h ai2 tt: h On h 612 

Ot ot ox ox 

= RuFiix, t) + RuGix, t) (23) 



and 



dFi(x,t) dG(x,t) , dFi(x,t) , (9G(a;,t) 
at ot ox ox 

= R2iFi{x, t) + i?22G(x, t), (24) 



where 



an = t, ai2 = 
021 = 0, 022 = t 



hii = —12AfX, 612 = ^^^^ (Sx — 3x In i 
621 = -l2AfX, 622 = -QAfPg{x) 



(25) 



(26) 
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and 

>. (27) 



R21 = 9Af (-X + fx In i) , R22 = 9A/ (ii - ^ + ln(l - x) + R'g{x] 

Eq. ()23j) and Eql24|l are two first order linear coupled differential equations. We now 
reduce these two equations to canonical form where they are decoupled from each 
other in terms of a new function. To that end we introduce a vector u{x, t) as 

G(M)) 

and write the two equations Eqs. ()2H|1 and ()24|1 in matrix form: 

aut{x,t) + bux{x,t) = Ru{x,t), (29) 
where the matrices a, b and R are given by 

(30) 
(31) 

\ "21 "22 / 

and 

R = [ 7 ) (32) 

y -0:21 K22 ) 

In Eq.(j2ni), Mt and are the derivatives with respect to t and x respectively. We 
note that the matrix a with its elements given by Eq. ()25|) is a non-singular matrix 
and hence multiplying Eq. (j29|l by from left, we get 

Ut{x,t) \ Aux{x,t) = B u{x^i^^ (33) 

where the new matrices A and B are: 

A = a-^b (34) 

and 

B = a-^R. (35) 

Equation ()33|) is a system of two coupled first order linear partial differential equa- 
tions in the two variables x and t for the vector u prescribed by Eq. ()28|) . Its principal 
part, i.e. Ut{x,t) + Aux{x,t) is completely characterized by the coefficient matrix 
A. Since the matrix A has n [here n=2] distinct eigenvalues, the system Eq.(|33l ) is 
a hyperbolic one and it is possible to obtain its canonical form in the following way 
HZIIIHI: 

Let Ai,A2 be the two distinct and real eigenvalues of the matrix A and let pi 
and P2 be the corresponding eigenvectors. Let P be a 2x2 matrix formed by the 
eigenvectors pi, p2, i.e. 

P=iPi P2). (36) 



an 


ai2 


^21 


0-22 ^ 


hi 


h2 


&21 


&22 


Rii 


Rl2 


R21 


R22 
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Now if A is the diagonal matrix with the eigenvalues Ai and A2 as the two elements, 
then we have 

P-'AP ^K=i\' M . 



A2 



(37) 



Let us define a new vector v by the relation 

v = p-^u (38) 

so that 

u = Pv. (39) 

For convenience, we have dropped the functional form in u{x, t), v{x, t) and P{x, t). 
Differentiating Eq. ljHTHl with respect to t and x respectively we get 

ut = Pvt + Ptv, = Pv^ + P^v. (40) 

Substituting Eqs.dSHI) and (001) in Eq.^, we obtain 

Pvt + Ptv + APv^ + AP^v = BPv. (41) 

Multiplying Eg. 1)411) from left by P^^ and using Eq. ()37)l . we obtain 

vt + Avcc = e, (42) 

where 

e = p-\BP-Pt-AP,)v (43) 



is a two component column matrix. Eq. ()42|) is in canonical form. In component 
form it is 

dvi ^ dvi . . „ / . . X 

2.2 Solution by the method of characteristics 

— , Ov ■ dv ' 

Equation (11^ shows that the principal part of the i th equation, viz. — ^ + Aj— 

ot ox 

involves only the component Vi of the vector v and its derivatives. To solve these 
equations by the method of characteristics, we define first the characteristic curve 
of the system .These are the curves in the x — t plane given hj x = Xi{t), where 
Xi{t) is a solution of the differential equation 

^ = Hx,t) (45) 

with Xi{x, t) being the eigenvalues of the coefficient matrix A{x, t) defined in Eq . ()34p . 
It is also known ^] that the characteristic curves of a hyperbolic system like Eq. ()33|) 
whose eigenvalues Xi{i = 1,2) are distinct, remain invariant under transformation of 
the system to its canonical form Eq. ()44|) . We also make the observation that along 
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a characteristic curve corresponding to an eigenvalue Aj, the left hand side ( i.e. the 
principal part) of Eq. lji^ is actually an ordinary derivative with respect to t, since 



along the characteristic curves. The actual integration on the right hand side of 
Eq. ljTTj) depends on the analytical solution of Eq. p^j) . It, in turn depends on the 
eigenvalues Ai and A2 of the coefficient matrix A{x,t) of Eq. ljHHj) . On calculating 
these eigenvalues and the eigenvectors with matrix a and b defined in Eqs ljHOj) and 
dnH), we find that they are too involved as discussed briefiy in Appendix. As a 
result, the characteristic equation cannot be solved analytically to find the 
curve X = Xi{t). 

The situation can be simplified by making some further approximations in Eqs. ()17l - 
EUj) . If we keep only leading terms a; In — compared to terms 0{x), then Eqs. ljlTH^UI) 

X 

become respectively 



dvijt) _ dvi{xi{t),t) ^ dvi{xi{t),t) dxjjt) 
dt dt dx dt 



(46) 





(48) 




(49) 




(50) 



and 



II{x,t)^R'g{x)G{x,t) + Pg{x) 



dG{x,t) 



(51) 



dx 



with 




(52) 




(53) 



The matrices b fEqlHTj) and R fEqlH^ are now given by 




(54) 



and 




(55) 
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The eigenvalues of the matrix A are then obtained from the characteristic equation 

det\A-\I\=0 (56) 

leading to 

Ai = (57) 

and 

A2 = ^. (58) 

The corresponding eigenvector matrix P fEqlHUjl is 

This simplification yields the column matrix e = ( ) defined in Eq. (j43|) and 



62 



occurred in Eq. ()44|1 with the components 



ei = \ ^- V2 + \ ^-^^ ^ + ^ ^« ^ j t;i 

and 



(60) 



_ -18A;n^x Ini /A^n^ 27A;n}a; In i 27A;n; / 11 1\\ 
^ 7t V ^ 49t 14t V 12 18 x/y ^' 

(61) 

Eqs.flHnjl and (jHH) will be used to solve Eq. ljTfj) for v\{t) and f2(i) on the character- 
istic curves corresponding to eigenvalues Ai and A2. 

Let (x, t) be a fixed point in the (x— t) plane through which the two characteristic 
curves X\{f) and X2(t) will pass i.e. 

Xi(f) = X and X2(t) = x. (62) 

Explicitly, the equations of the two characteristics fEqs. ()45|) ) are 

dxi{t) _ 9AfX Ini 
dt ~ t 

and 



(63) 



dx2{t) 2A/xlni 



(64) 



dt t 

satisfying the condition of passing through the point (x,f) as given by Eq. ()62|) . 
Solution of Eq. (j63|l is 

In Xi{t) = In X ( |j , (65) 
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while that of Eq.^ is 

In X2(t) = Inx Uj . (66) 

Furthermore, if the characteristic curves cut the initial curve t = to{= In ^) at 
xi = Ti and X2 = T2 respectively [Fig.l], then Eqs. (jU3j) and give 



and 



Inri = Inx ( (67) 



lnr2 = Ina; ( -^^ (6^ 



leading to 

r, = (69) 

and 

r^ = ^Cff''f\ (70) 
Dropping the bars over x and i , following are the expressions for the characteristics 

xi{t) = Ty^ (71) 

and 

X2{t) = T^'> . (72) 

Let us now find the solution of Eq.f lTfj) which has explicit forms on the characteristic 
curves as 

dvi{t) 



and 



ci {xi{t),t,vi{xi{t),t),V2ixi{t),t)) dt (73) 



= C2 {X2{t),t, Vi{x2{t),t),V2{x2{t),t)) dt, (74) 

Vi[t) 

where 

c,{x,{t),t,v,ix,it),t),V2{x,it),t)) = 12^/^1 (^) ^l f^^^T^l 

t Xi{t) \Vi{Xi{t),t) J 

( ISAfUf , 1 . / 11 nf , 1 M 1 

+ \ ^xi(t)ln—- + 9Af\ - + ln—-]}- 75 

\ 7 xi(t) V 12 18 xi{t)jjt ^ ' 

and 

18AfnfX2(t)ln^jTr 

C2 ix2it),t,Viix2it),t),V2ix2it),t)) = + {AfUf- 

27 , 2 / N, 1 27, /II nf\ 27, ,1 1 1 /^^i(a;2(t),t)\ 

— Am2.X2(t)/ra— - + —AfUf — + — AfUfln—- } - ) ) , i . 76) 

49 ^ / 2^ ^ X2 t 14 ^ ni2 187 14 ^ ^ X2 t J ^ U2 t ,^ / 
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Integration of Eqs. ()73|l and ()74|1 along the characteristic curves Xi{t) and X2{t) from 
to and i gives: 



to 



Ci {Xi (t) , t, Vl (Xi (t) ,t),V2 (xi (t) 



(77) 



V2{x,t) = V2{Ti)exp 



to 



Ci {X2 {t),t, Vl {X2 {t),t),V2 {X2 (t) ) dt 



{7t 



In Eqs. (l77p and (|78|) . we have now the boundary conditions Vi{ti) and V2{ti) at 
X = Ti and X = T2 where the two characteristic curves meet the initial line t = 
to(Figl). Eqs. lfTTj) and (fTSj) give the two components vi and t>2 of the function v 
defined in Eq. fl38|) . Now, from Eq. ljHHj) with found from Eq. ()59|) . we find that 
the components Vi and V2 are related to singlet structure function F2{x, t) and gluon 
momentum distribution G{x, t) by 



fi(x, t) = G{x, t), 



(79) 



and 



V2{x,t) 



x,t) - 



^G(x,t). 



(80) 



We can therefore rewrite Eqs. lffTjl and (fT8|l as 



G{x, t) = G{ti) exp 



to 



Ci {Xi {t') ,t',Vi {Xi (f) , t') , V2 {Xi {t') , t'))dt' 



^1) 



and 



3nf_ 
14 



G{x,t) + {Fi{T2)-'-^G{T2: 



exp 



/* C2 (X2 (f) , t', VI{X2 {t') , t') , V2 {X2 {t') , t') ) dt' 



12) 



In writing Eqs.(jH^ and (jH^ we have removed the bar over x and t and written the 
integration variable as t'. Eq. ()81|) and Eq. ()82j) are the most general solutions for the 
gluon and the singlet distribution within the present formalism. 

Unfortunately, the right hand side of both Eq. (j98|) and Eq. (j82p [see also Eqs. (j75|) 



and ()7(i|l ] contain the ratio y~j which is yet unknown. This ignorance forbids 

analytical forms of the quantities defined in both these equations. In order to get 
their analytical forms, we have to make additional plausible assumption about the 

ratio 



Vl 



. We assume that x and t dependence of vi and V2 are factorisable and 



the X dependent part does not deviate significantly from their values at the initial 
curve t = to, i.e. 

V2{x,t) V2{t2' 



vi{x,t) Vi{ti) 



■m, 



3) 
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where f(t) is some unknown function of t. Using Eq. (j83|l in Eqs. (j8T| and (j82|l . we 
get the gluon momentum distribution and the singlet structure function as : 



G{x, t) = G{Ti)exp 



2nf ( fioV-''^/ 



t 



1 



and 



X exp 



VIA 



^ V12 ■ 187 """to 



X \ \t ) 



9Af 



G[ti) J to xi{t'} 



Fi{x,t) = ^-^G{x,t)+{Fi{r,) 



14 



3iy 
14 



xexp 



G{t2) ] exp 



Fi{r2) - '-^G{t,) Jto 
where the function c'2{x2(t'),t') appearing on the r.h.s of Eq.(|H^ is given by 



9nf ix — x^ 



c'2ix2{t'),t')dt' 



15) 



c'2{x2{t'),t') 



27 
49 



27 



Afn,X2it)ln-^^ + --A,n, ^ 
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27 , , 

——AfUfln 



14 



^X2{t') t'f{t')- 



The functions xi{t) and X2{t) are defined in Eqs.([ZH) and ([Z2I). Eq.dHlI) and Eq.(jHSI) 
are the semi-analytical expressions for the gluon and the singlet structure functions 
in terms of the undetermined function f{t) . For analytical forms, we need to use 
simple test function for f{t) which we will discuss in ^ Further, if we expand f{t) 
about t = to and retain only the first term, then 

fit) = /(to) + (t - to)/'(t) +... ^ /(to) . (87) 

Integration over t' in Eqs. fl84|) and ()85|) can now be performed leading to the following 
analytical forms of gluon and singlet densities valid for t very close to the initial curve 
t = to in terms of one unknown parameter /o = /(to) : 



G{x,t) = Gin) exp 
X exp 



AFiin)-'-^Gin 



3 Gin) 
(^-(T''^')"4 



-./o + 



2nj 
7 



X 



exp 



U2 187 to 



and 



Fiix,t) 



^""f-Gix^t) + fe(r2) - ^G(r2) ) exp 



14 

X exp 
27 

Wf 1 x^ 

98 f 



GiT2) 



14 
1 



9Af[x- x^'-^^'^' 



n 



FiiT2) 



^G(r2) /. 



f. 27/11 
0/^^4^+14(12+1 



In- 



to 



X 



28 V V t 



2A, 



In— 

X 



19) 
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For phenomenological test we use our Eqs. (j84|l . (j85|) . (j88|l and (j89| . 



2.3 Compatibility with earlier solution 

In ^3] we derived an expression for the gluon momentum distribution considering 
only the gluon evolution equation and neglecting the quark singlet part. Therefore, 
putting F2 = in Eq. (jHH|) we should get back the earlier solution. But we find that 
the earlier solution cannot be exactly recovered from Eq. (jHH|) . The origin of this 
problem lies in the approximation Eq. ()53|) where we have taken Pg{x) ~ x Ini to 

be used in Eq.d^H). But in our earlier work we used Pg{x) as Pg{x) ~ x (\n ^ — j^)- 
Using that form of Pg{x) in Eq. ()54j] and without changing other quantities, the 
eigenvalues of the matrix A are found to be. 



Ai 



9AfX f-i| + lni 



(90) 



and 



A, 



-2AfX Ini 

t 



(91) 



to be compared with Eq. ljFTj) and Eq. (j^ . Only the first characteristic curve changes 
while the second one remains the same. 

Following the same procedure as in the derivation of Eq. (j88|) with the assumption 
Eq. ljHHj) and Eq.(jS7j) we find the the gluon distribution to be 



G{x,t) = G{T)exp[-H{x,t)] 
X exp 




^Ini 



(92) 



where the function H(x,t) is 
H{x,t) = 



12e 
- exp 



X 



11 
12 



In- 



G{r) 




(93) 



and r is given by 



exp 



, 1 11 
In - H 

X 12 



t 
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11 
12 



(94) 



which is exactly same as we derived in ^]. In Eq. ()93|) Fi and F2 are the incomplete 
Gamma functions 



Fi 



(95) 



12 



, 11 , 1 
0. + 



which is defined as 



r(o,z) 



(96) 



(97) 



Putting = in Eq.(|^ with the additional assumption that / (to) = 1, which 



f 2 

is equivalent to the assumption that the ratio — fEqITTj) along the characteristic is 



equal to its value on the initial curve, we find that Eq. ljU^ is reduced to 



G(x,t)=G(r) - 



1\ ^ t > fto\ 2 



t 



exp 



11 L_^_o}''^'' 



(98) 



which is exactly the same expression we obtained earlier in |^ without considering 
the singlet part. 



2.4 Non singlet structure function F2 {x,t) 

A similar procedure can be followed for the DGLAP equations for the non-singlet 
structure function ^2^"^ : 



di ~ T 
"1 dz 



{?> + Aln{l-x)}F^^{x,t) 



Following the approximation similar to Eqs. (jl5|) and (jl6|) we get, 



F^\-,t)^F^'{x,t)+[xY. 



u 



k=l 



dF^^{x,t) 
dx 



(99) 



(100) 



Now putting Eq.JTUni) in Eq-dHHI) we get 

5^2^^ (x,t) A 



dt 



{3 + 4/n(l - x)} Ff^ix, t) + 2 -^{z^ - l)Fi'^x, t) 

Jx 1 — Z 



+2 



1 dz 



: l-Z 



;i+^')(^E 



u 



k=l 



dFf'{x,t) 
dx 



.(101) 



Carrying out the integration in z and neglecting terms O(x^) and higher, we can 
write Ea. (ITTTT|l as 

dF^^{x, t) SAf X dF^^{x, t) Af {4/n(l - x) + 2x} 



dt 



3 t dx 



F^''{x,t) (102) 



Eq. p02|) is a partial differential equation for the non-singlet function F,^^{x, t) with 
respect to the variables x and t. The characteristic curve of the Eq.( I102|l is given 
by the solution of the differential equation 



dx^{t) 



^Afx 



dt 3 t 

13 



(103) 



Assuming that the characteristic curve passes through a point (x, t) i.e X3{t) = x in 
the X — t space, we get the solution of Eg . (11031) to be 

X 3 t ^ ' 

If the characteristic curve cuts the initial curve t = to at a, point x^lto) = t^, then 
Eq. (|TTH|l gives 

ln^ = -^lnh, (105) 



so that 



r3 = 5(T) • (106) 



t 

Dropping the 'hats' over x and t, the equation of the characteristic is 

to\ — 



xsit) = r3[j) . (107) 

On using Eq. p03|) in Eq. ()102|) , the left hand side becomes an ordinary derivative 
with respect to t and the equation becomes an ordinary differential equation: 

^^^""'^f'^''^ = c^Hx{t),t)Fi^Hx{t),t) , (108) 

where 

Integrating Eq. ()108|) over t from to to i along the characteristic curve x^it) fEq llOTj) 
and finally dropping the bars over x and t we get the solution for the non-singlet as 



X,t) 




(110) 



Eqs. (l88p . (|89|) and (jliop are the analytical solutions of the DGLAP equations within 
the present formalism. While Eqs. (j88j) and (j89|) contain a free parameter /q, Eq. (jllO|) 
is a parameter free solution. 

Using our results derived in this section, we will calculate the proton structure 
function Fj' which is related to the singlet and the non-singlet structure function by 
the equation 

= ^Fi'' + -Fi . (Ill) 
^ 18 ^ 18 ^ V ; 

We discuss in the next section the phenomenological consequences of our results 
derived in this section. 
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10"^ 



X 

Figure 1: Characteristic curves given by Eq. (f7T|) fsohd hne),Eq. ((721) (dotted hne) and 
Eq. ((TOTj) (dash-dotted hne) drawn from the initial hne t = to{= luj^) to the point 
X = 0.001, i = 7.65(Q2 = SOGeV^). Here we take Uf = 4:,Af = ^, Ql = IGeV^ and 
A2 = 0.200GeV'^. 

3 Results and discussion 

In §2.2l we have derived semi-analytical [Eq. (|H^ and Eq. ljHHjl ] and analytical [Eq. (|HH|) 
and Eq.(|HU|)] solutions of the gluon and the singlet and in ^2.4l obtained an analytical 
solution [Eq. ()110|) ] for non-singlet structure function ^2^"^ in LO.Let us now compare 
them with the exact solutions of the DGLAP equations to see in what region of x—Q^ 
space they agree with the exact ones within a few (say 15%) percent. We choose 
MRST2001 LO ^ solutions for comparison. For evolution of Eq.((HlD, Ea. (jHH|l . 
Eq. ljHHjl . Eq. ljH^ and Eq. (jllOj) we take the inputs from the above refdH]. We note 
that in our expressions [Eq.dHl), Eq.(IH3), Eq.dHHI), Eq.dHHI) and Eq.jnU))] we need to 
take the inputs as functions of Tj (i=l,2,3) rather than of x, where are the points 
at which the characteristic curves defined in Eqs. (j71|) . (j72j) and (j98p cut the initial 
line t = to {t = In ^) respectively. For nj = A and Aj = (^/ = si")' ^^ese are 

36 ^ 32 

Ti = x^~) " ■, T2 = x^~) and = x^~^ ^ respectively. The inputs thus acquire a 
t dependence also in addition to x and are obtained by a formal replacement x ^ Ti 
in the MRST2001 LO ^[9^ input distributions. 

To evaluate numerically the semi-analytical expressions Eq. ()84|) and Eq.()85p. we 
need the unknown function f(t) defined in Eq. (jSH|) . As f(t) is the ratio of the t 
dependent part of V2 and Vi defined in Eqs. lfTUI) and (jSHI) we assume it to be of the 
form fit) = (Int)'^, where is a parameter to be determined by comparison with 
the exact solutions. Varying k continuously from negative to positive values, we 
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Figure 2: Semi-analytical singlet structure function F2{x,t) fEq lH^j) compared with 
the exact MRST01LO^19j solutions (solid line) at nine different for three different 
values of the parameter k = 1.0 (dashed line ), 1.2 (dotted line ) and 1.4 (dash-dotted 
line.) 
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Figure 3: Semi-analytical Gluon momentum distribution given by Eq. ()84|) compared 
with the exact MRST01LO|Tn] solutions at nine different representative Q^. The 
gluon (EqlHH) is shown for only one value of = 1.2 as discussed in the text. 
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find that only in a limited range of positive values from A; ~ 1.0 to /c ~ 1.4, the 
semi-analytical solutions G{x,t) [Eq. ljH^ ] and F2{x,t) [Eq. (jH3|) ] compare well with 
the exact ones. In Fig. 2 we show F2{x, t) [Eq. (j85|) ] as function of x at nine different 
representative for three different values oi k = 1.0, 1.2 and 1.4. We see that the 
value of k is crucial for the consistency of our solution vis-a-vis the exact ones and 
different values of k seem to suit at different Q^. Our prediction remains below the 
exact solution for up to ^ 20GeV^ and x < IQ-^. But above > QQGeV^ 
and X < 10~^, it overshoots the exact solution and rises faster with decreasing x. 
However, above x > 10~^, our prediction lies above the exact ones for all the 
values shown. Taking the graph for k = 1.2, we calculate the percentage deviation of 
our solution from the exact one at seven different equally spaced x at each and 
find that in the range 6.5 GeV"^ <Q^< GOGeV"^ and 10~^ < a; < 10~^, our solution 
(Eq.(jHSI)) is within 15% of the exact one. 

Taking = 1.2 we evolve the gluon given by Eq. ()84|) and compare with the exact 
MRSTOl ^H] solution in Fig. 3 at nine different Q^. We see almost similar results 
as in the case of ^2^. Below ~ 20 GeV^ the predicted gluon density remains 
lower than the exact one in the small x {x < 10~^) region but as > 20GeV^, it 
begins to overshoot the exact solution, so that at ~ 60GeV^ and x ~ 10"'^, the 
percentage deviation from the exact one becomes more than 15%. For x > 10~^ 
and 5 GeV"^ < < 60 GeV"^, the predicted gluon(EqlHH) is within 15% of the exact 



In Fig. 4 the non-singlet solution F^^{x,t) fEq. ()110|) ) is shown along with the 
exact MRST2001 LO solution. Unlike the singlet and the gluon, the non-singlet 
solution is free from any free parameter. The qualitative features of Eq. (jllOj) and 
the exact solution are identical, the exact one remaining always slightly higher. 

Combining Eq. ljH^ and Eq. ()110p in Eq. ()lll|) we calculate the proton structure 
function . The result is compared with the exact solution jTHI in Fig. 5 at twelve 
different representative values of as a function of x. Taking some (eight) equally 
spaced x values at each Q^, we calculate the percentage deviation of our result from 
the exact one and find that for 5GeV'^ < < AbGeV"^ and 10~^ <x< 10~^, our 
F| is within 15% of the exact MRSTOILO solutions. 

Next we study the analytical solutions Eq. ()88|) and Eq. ()89|) which contain a free 



parameter /q and are expected to be valid at t very close to to ( = 1^ tt^ ) • Here also 



we find the value of the free parameter /o such that our solutions are compatible 
with the exact ones in a definite x — Q"^ range. We find that for /o = 0.8, our 
approximate analytical solutions reproduce the exact results in a limited range of 
a;(10-^ <x< 10-2) and Q'^{hGeV'^ < < QOGeV^) within 15%. 

4 Conclusion 

We have presented in this paper solutions of the singlet structure function F^i^x^t) 
and the gluon momentum distribution G{x,t) valid to be at small x. Applying the 
method of characteristics , the LO coupled DGLAP evolution equations are solved 
in the x — t space. The results are presented both in analytic and semi-analytic 



solution. 
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10"' 10"* 10"' 10"' 10"' 10"' 10'* 10"' 10'" 10"' 10"' 10'* 10"' 10"' 10"' 



X X X 

Figure 4: Non-singlet structure function F.^^ (Eq llin|l compared with the exact 
solution at the same values of as in Fig.Q and Fig.©. 
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Figure 5: Proton structure function F^x, t) calculated using Eq. (jHH|) and Eq. ljllOj) 
in Eq. ljllH) compared with the exact MRSTOILO PHj solution. Only one value of the 
free parameter k{= 1.2) is used as discussed in the text. 
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forms and compared with the exact MRST2001LO solutions to find the range 
of vahdity of the solution. Applying the same method, LO evolution equation for 
the non-singlet structure function F:^^{x,t) is also solved analytically. 

In order to obtain the solutions for F2 and G, we however need additional infor- 

mation about the ratio — = — ^ ^.This is achieved through an unknown 

vi G{x,t) 14 

function f{t) fEq l83|) or an unknown parameter fn fEq lHTj) . Through their suitable 
choices, f{t) = {IntY''^ or /o = 0.8, we found regions oi x — where our approxi- 
mate solutions are compatible with the exact ones within 15 %. While this additional 
information needed is an inherent limitation of the present formalism, still it can be 
considered an improvement over our previous work in refs |il3ti22t i23i| where identical 
t— evolution were assumed for gluon and singlet distributions unlike Eq. (jH^ . 

In this paper we have not compared our results with data directly, since exact 
LO solutions like JH] reproduce the data correctly.A similar agreement with data 
in the entire x — Q^ range of HERA of our approximate solution presumably need a 
new fit of the input distributions and other free parameters/functions like f{t) and 
/q. The same inputs as used [19^ by us to compare our results with exact ones are 
not sufficient to reproduce the entire data, as is evident from our analysis. This is 
understandable, since approximate solutions refer to a different procedure with dif- 
ferent evolution and thus to a possibly different input. Such a possibility is currently 
under investigation. 
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